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l-l ^Abstract 

' ' .The core mass of Saturn is commonly assumed to be 10-25 M^ as predicted by interior models with various equations 
I^ of state (EOSs) and the Voyager gravity data, and hence larger than that of Jupiter (0-10 M^). We here re-analyze 

Saturn's internal structure and evolution by using more recent gravity data from the Cassini mission and different 

Ojphysical equations of state: the ab initio LM-REOS which is rather soft in Saturn's outer regions but stiff at high 

^ pressures, the standard Sesame-EOS which shows the opposite behavior, and the commonly used SCvH-i EOS. For all 

three EOS we find similar core mass ranges, i.e. of 0-20 M® for SCvH-i and Sesame EOS and of 0-17 M^ for LM-REOS. 
^ Assuming an atmospheric helium mass abundance of 18%, we find maximum atmospheric metallicities, .^atm of 7x solar 

for SCvH-i and Sesame-based models and a total mass of heavy elements, Mz of 25-3OM0. Some models are Jupiter- like. 
I— lAVith LM-REOS, we find Mz = 16-2OM0, less than for Jupiter, and Z^tm ^ 3x solar. For Saturn, we compute moment 
Ph of inertia values A = 0.2355(5). Furthermore, we confirm that homogeneous evolution leads to cooling times of only ^ 2.5 
W iGyr, independent on the applied EOS. Our results demonstrate the need for accurately measured atmospheric helium 
r^ and oxygen abundances, and of the moment of inertia for a better understanding of Saturn's structure and evolution. 
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^ 1. Introduction 

. Saturn is the planet with the lowest mean density in 
►^ "the solar system. Since the mechanisms that can in- 
r — flate exoplanets with observed overlarge radii do not hold 
f^ for the outer planet Saturn, one might thus intuitively 
C .think of Saturn as having a smaller core and smaller over- 
'all metallicity than Jupiter. However, quantitative esti- 
mates on the core mass and on the total heavy element 
fenrichment solely come from interior model calculations, 
and the same modeling approach applied to both plan- 
. . fets just predicts the opposite: an about two times larger 
_ ^ maximu m core mass and he a vy el e ment enrichment for 
K^ Saturn (jSaumon and Guillotl . l2004l ICuillot and Gautier . 
^ I2OO7I ). A higher envelope metallicity of Saturn is also sup- 
.. ."ported by the measure d atmospheric C:H ra tios, which is 



o 



2009; scaled to the 



20031 ) but only 



'^ 9x solar for Saturn ([Fletcher et al . . _ 
solar system abundanc e data of LoddersI 
3-5 X solar for Jupiter (JAtreva et al.l . 120031 ). 

Certainty about the present core mass and enve- 
lope metallicity is desirable because these parameters 
contain information — albeit not ne cessarily uniquely 
(JHelled et al.l I2OIOI: iBolev et al.l . l201l[) — on the forma- 
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tion environment, i.e. on the protosolar disk, and on the 
process of formation. 

Models bv lSaumon and Guillotl (|2004l ). hereafter SG04, 
are often considered the standard of what we know today 
about Saturn's present internal structure i n terms of core 
mass and heavy element enrichmen t (e.g., lAlibert et al.l . 
2OO5I : IPodson-Robinson et al.l . I2OIOI) , for mainly two rea- 
sons. First, these models have been computed for various 
physical equations of state(EOS) for Saturn's likely main 
constituents H and He that also give acceptable solution 
for Jupiter's interior and evolution (the EOSs SCvH-i, LM- 
H4, LM-SOCP). Independent on the EOS, the possible 
core mass range was found to be ~ 1O-25M0, while for 
Jupiter ~ O-IOM^. Second, a wide range of input param- 
eters was accounted for such as a the position of an internal 
layer boundary that separates a helium-poor, outer from 
a helium-rich, inner envelope. However, SG04 computed 
constant metallicity envelope models only, an assumption 
that tremendously restricts the resulting range of interior 
models. 

In earlier m odels by Gudkova & Zharkov (1999) and 
Guillotl (|l999l) . the metallicity was allowed to vary across 



the internal layer boundary. As a consequence, zero-core 
mass models with high heavy element enrichment in the 
deep envelope were found for both Jupiter and Saturn. 

The new Cassini gravity data with their tight ob- 
servational error bars, and also l ong-term observ a tional 
data of the Saturnian system ([Jacobson et al.l . 12006 : 
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Anderson and Schubertl . 120071 ) raised hope to better con- 
strain Saturn's internal structure. Surprisingly, the most 
recent Saturn models based on those gravity data cover 
an eve n bigger, minimum core mass range of '^ 0- 



30Mm (lAnderson and Schubertl. [2007l:lHelledet al.l.l2009at 



Helledl . l201l[) . Therefore. iHelledl (|201l[ ) suggests to mea 



sure the axial moment of inertia as an additional con- 
straint. Her models, however, employ empirical pressure- 
density relations that may reach out of the realm of physi- 
cal EOS which agree with the ava ilable experimental data 
(see, e.g. SG04: iHolst et al.ll2012h . 

Our Saturn models are the first that are based on both 
physical equations of state and the Cassini data. Not is 
it the purpose of this work to better constrain the core 
mass: this cannot be achieved within the standard three- 
layer modeling approach, which is adopted in this work. 
Instead, we here investigate the overall behavior of core 
mass, atmospheric metallicity, and deep envelope mctal- 
licity on the input parameters: we vary the position of 
an internal layer boundar y in order to recall its influ ence 
on the core mass, see also lCuillot and Gautien ( 2007 ): we 
exchange the EOS of the envelope material (LM-REOS, 
SCvH-i EOS, Sesame EOS), and we adopt two different 
periods of rotation of lOh 32m and lOh 39m. In lack of 
accurate observations, we make predictions on the possi- 
ble helium and heavy element mass fractions in Saturn's 
atmosphere in dependence on the J4 value and the un- 
certainty in the rotational period. Our results on the at- 
mospheric helium abundance can serve as constraints for 
future models of He-sedimcntation in Saturn, as long as 
Saturn's atmospheric IIe:Il2 ratio is not accurately mea- 
sured. 

Observations of young stellar systems and protostellar 
disks commonly point to forrnation of the giant planets 
within a few Myr (jStrom et al.l . 119931 ) , implying a billions- 
of-years-old planet should have the same age as its host 
star. However, homogeneous evolution calculations for 
Saturn, which are mainly based on t he SCvH-i EOS, gen- 



erally yield co o ling t imes of 2-3 Gyr (jSaumon et al.l . 11992 



Fortnev et al.l . 1201 ll ). about only half of the age of the 



Sun. This implies a higher luminosity of present Saturn 
than it should have if the underlying assumption of homo- 
geneous evolution would hold. Despite the obvious failure 
of this assumption, we here adopt it once more in order to 
investigate the influence of the EOS on the cooling time. 
In Section F2.1l we describe our modeling procedure. Sec- 
tion 12.21 is devoted to a detailed description of the obser- 
vational data, and Section 1^^ to the applied EOSs. Our 
results are presented in Section [3] In Section 13.11 we in- 
vestigate the influence on different H-Hc-EOS on Saturn's 
structure and in Section 13.21 of the atmospheric He abun- 
dance and rotation rate. In Section 13.31 we give the val- 
ues for the non-dimensional moment of inertia. Section 
13.41 contains the cooling curves. Section Hj includes a dis- 
cussion on the implications for Saturn's formation process 
(|4.3p . on the applicability of the three- layer assumption in 
the presence of He rain (|4.4p , and a summary of our main 



findings (|4.6 



2. Methods 

2.1. Planetary structure modeling 

For understanding the interior of giant gas planets like 
Saturn it is necessary to consider the gravitational field of 
the planet. The shape of the field is influenced by differ- 
ent effects. Saturn for instance has primarily the form of 
an ellipsoid due to its rapid rotation, which can be seen 
from the rather high ratio of centrifugal to gravitational 
forces, q = u'^ R^J{GM), where uj is the angular veloc- 
ity, i?eq is the equatorial radius, and M the total mass. 
For Saturn, q ^ 0.155 with an uncertainty of 0.004 due 
to the uncertainty in the rotation period and equatorial 
radius (see Section \T^ . for Jupiter, q = 0.089, and for 
the Sun, q = 0.00002. Tidal forces caused by the gravity 
of the moons or the parent star can also change the form 
of a planet's gravity field. While this effect can be impor- 
tant for close-in exoplanets it is tiny for Saturn and has 
not been measured yet for any giant planet in the solar 
system. To assess the rotationally induced deformation, 
the gravity field $''^'' exterior to the mass M is expanded 
into a series of Legendre polynomials P2n, where the ex- 
pansion coefficients J2n are the gravitational moments at 
the equatorial reference radius i?oq, 



J2 



1 



MRl^ 



d:\p{r,eyP2n{t) 



(1) 



Being integrals of the internal mass distribution over the 
volume enclosed within the geoid of equatorial radius iieq, 
the J2n can be written as depth-dependent functions J2n(0 
whose values increase continuously from the center out- 
ward until the observed values Jj" are reached at the 
geoid's mean radius I ~ i?„i. As a measure for the contri- 
bution dJ2n of a shell at / and extension dl to j!^^ ^ we 
can define the normalized contribution function 



C2n(0 



{dJ2n/dl)\l 
Jdl' {dJ2n/dl') 



(2) 



For modeling Saturn w e use the same method and code 



Nettelmann et al.l (J2012r ) for Jupiter. We adopt the 



standard three-layer structure with two envelopes and a 
core. The composition of each of the envelopes is diverted 
into the three components hydrogen, helium, and heavy el- 
ements, whereas the core consists of heavy elements only. 
The helium mass fractions and the metallicities (i.e. the 
heavy element mass fractions) are parameterized by Yi, 
Zi and Y2, Z2 for the outer and the inner envelope, re- 
spectively. This implies the assumption of homogeneous 
envelopes. The transition between them occurs at the 
transition pressure Pi_2 which is a free parameter. As ob- 
servational constraints we take into account i?cq, w, the 
total mass Msat, the temperature Ti at the 1 bar level of 
the planet, and the lowest order moments J2 and J4. 



For given values of Yi and of the mean helium abundance 
Y, I2 is adjusted to fit Y, while Zi and Z2 are adjusted 
to fit J2 and J4 . Mass conservation is then ensured by the 
choice of the core mass Mcorc = "^(^core)- 

2.2. Observational constraints 

While the Cassini mission could provide tight con- 
straints on Saturn's gravity field, there are still important 
remaining uncertainties, in particular in Saturn's period 
of rotation, equatorial radius, and the atmospheric helium 
abundance. 

Period of rotation. Prior to the Cassini observations, Sat- 
urn's period of rotation was taken to be lOh 39m 24s, 
the detected periodicity in the kilometric radio emissions 
of Saturn's magn etic field as measured by the Voyager I 
and II spacecraft ( Desch and Kaiseil Il981l) . Cassini how- 



a reference geoid at 1 bar to be ^ 1 00 km smaller than 



ever revealed a prolongation of this period by several min- 
utes within just 20 years; thus the observed magnetic 
field modulati ons may not re f lect t he rotation of Saturn's 



deep interior (jGurnett et al.l . 120071) . On the other hand 



while alternative methods of deriving the rotation rate 
from observed wind speeds make assumptions that may 
not hold true, such as the minimum energy of the zonal 
winds or a minimum heig ht of isobar-surfaces rela t ive to 
computed geoid sur faces (lAnderson and Schubert! 12007 : 



Helled et al.l . l2009bl ). that alternative methods just sug- 
gest similar values of ~ lOh 32m. We therefore use these 
values as the uncertainty in Saturn's real solid body rota- 
tion period and compute interior models for both periods, 
i.e. for lOh 32m and lOh 39m. Note that we neglect here 
the uncertainty to Saturn's structure from the possibility 
of differential rotation on cylinders. On the other hand, all 
observational wind data can well be reproduced by th e as- 
sumption of solid-body rotation (JHelled et al.l . l2009bll and 
the effect of zonal winds on Jupiter's structure have been 
shown to be negligible if thei r penetration depth is limited 
to 1000 km (|Hubbardl . ll999l) . 



Equatorial radius. Our method of interior modeling re- 
quires all outer boundary conditions to be provided at the 
same pressure level. Although the Voyager I, II, and Pi- 
oneer 11 observational data refer to the isobar-surface at 



100 mbar (jLindal et al.l . 119851) . we prefer to use the 1 bar 



level as outer boundary. In particular, we use the equa- 
torial 1 bar radius of Rr. 



the given value in lLindal et al.l (|1985I ). Neglecting the de- 
scribed inconsistency, we here use this i?eq value, and we 
use it independently on the period of rotation value. To 
our awareness, the only consistent reference systems for 
both limiting periods of ro tation ( 1 Oh 3 9m 24s and lOh 
32m 35s) are presented in iHelledl ( 20111 ). Her reference 
system values arc listed in Table [T] for comparison. The 
surface radius -Rsat is the mean radius of the reference 
geoid with the equatorial radius i?oq- 

Gravitational coefficients. In order to get the gravity 
field (J2, J4, Je) at the outer boundary as defined 
by the equato rial 1 bar radiu s of 6 0,268 km, we scale 
the values of Ijacobson et al.l (J2006r) . which arc Jj = 
16290.71(0.27), 4 = -935.8(2.8), and J^ = 86.1(9.6). 



We use the scaling relation Ji^^Rcq 



J2nRc 



where 



R, 



' 2n 



Rq 



60, 330 km is the equatorial reference radius in 



Jacobson et al. 



Jacobson et al. 



( 2006 ). The gravitational coefficients of 
(I2OO6I ) arc based on, but not limited to, 
Pioneerll, Voyager, and Cassini tracking data as well as 
long-term Earth-based and HST astrometry data. They 
have significan tly reduced error bars compared to the Voy 



19891 ). Si nce the 



ager era data (iCampbell and Anderson 

- Reg v alue used in this work is larger than that in iHelled 

(|201l[) . the scaled, absolute J2„ values are smaller. Under 



1985t ICuillot and Gautierl . 120071) . This radius is a com- 



variation of the J4 value within its 6 a erro r bars ( Sectio n 
13. 2p . we also cover the J4 values used in iHelledl (|201ll ). 
see Table HJ Therefore, the obtained sets of solutions can 
reasonably be compared to each other. 

Mean helium abundance. For the protosolar cloud where 
the Sun and the giant planets formed of, Bahcall & Pin- 
sonneault (1995) calculate a mean helium abundance of 
27.0 to 27.8% by mass, depending mainly on the inclusion 
of helium and heavy element diffusion into solar evolution 
models. We require our models to have and mean helium 
abundance Y = 0.2750(1) by mass with respect to the 
H/He subsystem. By particle numbers, this corresponds 
to He:H2 := A^hc/^h^ = 0.1732. 

Atmospheric helium abundance. Combined Voyager in- 
frared spectrometer (IRIS) data and Voyager radio occul- 
tation (RSS) temperature profile data revealed a deple- 
tion of He in the atmospheres of b oth Jupiter and Saturn 



60,268 km ([Lindal et al.l . compared to the protosolar value (jConrath et al.l . Il984[ ). 



puted one of a geoid that rotates with a solid body period 
(System HI period) of lOli 39m 24s plus an additional, lat- 
itude ((/))-dependent component according to the observed 
zonal wind speeds. The difference in radius at the equa- 
tor to a reference geoid for a rigidly rotating Saturn in 
hydrostatic equilibrium — the appropriate radius to con- 
strain interior models — is ~ 120 km at the 100 mbar level. 
The exact di f ferenc e at the 1 bar level is not provided (see 



Lindal et al.l . Il985l for details) , but we expect it to be of 



same size, and therefore would expect the R^q value of 



In particular, a modest depletion He:H2 = 0.110 ± 0.032 
was found in Jupiter, and a strong depletion He:H2 = 
0.034 ± 0.024 in Saturn. These particle ratios correspond 
to atmospheric mass mixing ratios yatm.j = 0.18 ± 0.04 
for Jupiter and Yatm, s = 0-06 ± 0.05 for Saturn. How- 
ever, the He abundance detector (HAD) aboard the Galileo 
probe measured in situ Jupiter's atm ospheric He:H2 to 
be 0.157 ± 0.003 (|von Zahn et al.l . ll998l) . corresponding to 
^atm, J = 0.238 ± 0.006. Because of the discrepancy to the 
former results, the Voyager data for Jupit e r and Saturn 
were re-examined by IConrath and Gautieii (|2000| ) . They 



found the Voyager data for Jupiter could be made con- 
sistent with the Gahleo data if the temperature profile 
obtained by Voyager RSS measurements were shifted by 
2 K towards colder temperatures. Because of a possible 
system atic error of the RSS data, IConrath and Gautiei 
( 20001 ) developed an inversion algorithm to infer the He:H 
ratio and the temperature profile from the Saturnian IRIS 
spectra alone. Their results suggest a significantly larger 
He abundance yatm,s = 0.18-0.25, in agreement with the 
early Pioneer IR data based determination of 0.18 2 ±0.005 
(He:H2 = 0.111 ±3%) (lOrton and Ingersolill98n|). U nfor- 



tunatcly, the method of IConrath and Gautieii (J200 



not be checked by application to Jupiter due to disturbing 
NH3 cloud formation in the spectral range of interest in 
J upiter's s hghtly warmer atmosphere. 

KerlevI ( 2004al) instead suggests to trust the ratio of 



the He abundances for the two planets as derived orig- 
inally from the Voyager IRIS and RSS data. With 
(Hc:H2)s = 0.31(1 ± 1.1) x (HciHajj and mne = 4 g/mol, 
TOH2 = 2 g/mol, this would give 

_ (Hc:H2)s mH2 +mHe(He:H2)j 
(He:H2)j mH2+mHo(He:H2)s 

which is in numbers 



Y^tm, s = 0.353 (1 ± 2) X 0.238 = 0.084 (1 ± 2) 



(4) 



We thus compute Saturn models for different helium abun- 
dances Yatm between 0.10 and 0.18. Because of the as- 
sumed convection below the 1 bar level, Yi = Yatm- 

Atmospheric temperature. Saturn's atmospheric tempera- 
ture was determined to be 135 ib 4 K at the 1 bar level 
and 145 ± 4 K at 1.3 bars (JLindal et all . Il985[ ). For our 
model calculations we assume a slightly higher 1 bar tem- 
perature Tibar of 140 K but also compute single models 
with Tibar = 135 K. Note that the physical tempera- 
ture is not a direct observable. In particular, the obser- 
vational constraint T i bar d epends on an assumed com- 
position (JLindal et al.L Il985r) . It is the height-dependent 
refractivity which was measured during egress and ingress 
of the Voyager II spacecraft. From these data and the as- 
sumed composition, the mass density and thus the particle 
number density can be derived. Integration of the equa- 
tion of hydrostatic equilibrium, dP/dh = —g p, over height 
h in the atmosphere allows to relate density to pressure. 
Finally, the thermal equation of state for the assumed com- 
position with mean molecular weight p, yields the temper- 
ature T{P,p), which is proportional to JlP/p. Since the 
abov e cited Ti bar value was derived for the low Y value of 
0.06 (JLindal et al.L 119851 ). a revision of that value may also 



require a revision of the atmospheric temperature determi- 
nation towards higher temperatures. Table [T] summarizes 
the used constraints. 



Tabic 1: Observational constraints for Saturn 



Parameter 



default value 



Helled (2011) 



P = 27r/w 


lOh 39m 
lOh 32m 


lOh 39m 24s", 


lOh 32m 35s'= 


Afsat (lO^^g) 


5.683566'^ 


- 


- 


Req (km) 


60268^ 


60141.4 


60256.9 


■h (10-4) 


16324.2(0.3) 


16393.1 


16330.2 


J4 (10-4) 


-939.6(2.8) 


-947.6 


-940.4 


Je (10-4) 


86.6(9.6) 


87.8 


86.8 


Tlbar (K) 


140 K 


- 


- 


Teff (K) 


95.0 ±0.5^= 


- 


- 



iDesch and Kaisej | |1981| ): from kilometric radiation and magnetic 

fi eld data 

Lindal et aj] lll985l'li for P=10 h 39m 24s 



Anderson and Schubert] d2007l ): lOOmbar isosurface height minimiza- 



tion 

'^ with 



GMsat from iJacobson et al.l | |2006|) and G 
6 .67384 IQ-i'^m^s-^ (COD ATA 2010) 
" iGuillot and Gautied | |2007|) 



of the set of possible structure models for present Saturn, 
wc pick one model, implying known values of the structure 
parameters Mcoro, Y2, Zi, and Z2. Keeping these param- 
eter values constant, we then generate a series of models 
with increased Ti bar values. For these models, we also 
keep the angular momentum, L conserved and store the 
energy of rotation, E-^ot = Lu/2. The higher Tibar , the 
warmer the interior and hence the larger the planet ra- 
dius Rp and the higher the luminosity. Finally, Tibar is 
mapped onto time by integrating the cooling equation 



f-M 



dt = - 



/^r"* dm T[t) ds[t) + M,orcCvdT,orc{t) + dE,,t{t) 



4nR^{t)<j{T^\{t) 



mt)) 



^radioi,t;j 

(5) 
from present time tq = 4.56 Gyr backwards. In Equa- 
tion ([5]), parameters that change with time are written 
as a function of time, although most of them do not de- 
pend on time explicitly; in fact, only the solar luminos- 
ity Lq, which is proportional to T^ and assumed to in- 
crease linearly with time, and the luminosity Lcorc from 
the d ecay of radioactive eleme nts in the core do so. Un- 



like in lNettelmann et al.l ()2012l) where Jupiter's iradio was 



kept constant, we here include its time-dependence as in 



Nettelmann et al.l (|201ll ). because Saturn's Mcoic/M ratio 
can be ~ 6 times larger than Jupiter's and thus special 
core contributions potentially be non-negligi ble. The rela- 



tion b etween Toff and Tibar is adapted from IGuillot et al 
(|l995l) with constant K ~ 1.565. As Tibar determines Rp 
and Tcfi, the map Tibar '-^ t includes the more familiar 
maps i?p i-> i and Tcff i~> t. We then repeat the described 
procedure to compute the cooling curve for another struc- 
ture model of present Saturn in order to learn about, e.g., 
the effect of the chosen EOS on the cooling time. 



2.3. Planetary evolution modeling 
We compute the cooling t ime of Saturn almost exactly as 



Nettelmann et al.l ()2012l ) for Jupiter. In particular, out 



2.4. Equations of state 

Saturn's mantle is believed to mainly consist of hydro- 
gen and helium, and some amount of heavier atoms or 



molecules which we call heavy elements. For each of these 
three components we apply a separate EOS. To obtain an 
EOS for the mantle material, we linearly mix the single- 
component EOS. Saturn models are then computed for 
different sets o f equations of state for Saturn 's mantle: 



LM-REOS (see iNettelmann et al.l . l2008l l2012l and refer 
ences therein) SCvH EOS for hydrogen and helium by 



Saumon et al.l ( 19951 ) where we mimic heavy elements by 
scaling the density of the He EOS by a factor of 3/2, and 
the Sesame EOS. 



Th e Sesame-5251 hydrogen EOS (|Lvon and Johnsonl . 
19921 ) is the deuterium EOS 5263 scaled in density as de- 
veloped by Kerley in 1972. It is based on the chemical 
picture and built upon the assumption of three phases: a 
molecular solid phase, an atomic solid phase, and a fluid 
phase which takes into account chemical equilibrium be- 
tween molecules and atoms and ionization equilibrium be- 
tween at oms, protons and electrons. A completely revised 
version (IKerlevl 120031 ) includes fits to more recent shock 
compression data resulting into a larger compressibility at 
^ 0.5 Mbar, and a smaller one at ^ 10 Mbar. We here 
apply the earlier version as it shows stronger deviations 
to the SCvH EOS and the LM-R EOS, which allows to 
attribute differences in the resulting Saturn models more 
clearly to properties of the EOS. For application to plane- 
tary models, th e Sesam e H EO S is linearly mixed with the 
helium EOS of lKerlevI (|2004b[ ) and with the water EOS 
H2O-REOS. 

2.4- 1- Core material 

The cores of giant planets are usually assumed to con- 
sist of ic es and rocks a s a relict from their time of for- 
mation (JMizunol . Il980t iMiller and FortnevL I2OIII) . As 
limiting cases, we either compu te models with pu r e wa- 



ter cores using the H2O-REOS (jNettelmann et all . 12008 



French et al.l . l2009l). or with pure rocky core s using the 
P — p relation of iHubbard and MarlevI ( 19891) . This ap- 
proach ignores the possibility of an eroded core that would 
also contain some hydrogen and helium. 

2.4-2. Comparison of the applied hydrogen EOSs 

Even if the heavy elements in Saturn would have the 
low molecular weight of helium, still at least 70% of the 
particles in Saturn would be hydrogen. Therefore, among 
the equations of state for planetary materials, that of hy- 
drogen is expected to have the biggest influence on the 
resulting structure models. In Figure [T] we compare the 
above described hydrogen EOSs with experimental deu- 
terium shock compression data. For a certain experimen- 
tal set-up and given initial conditions (/9o,2o), different fi- 
nal compression ratios p/ po and pressures can be achieved 
depending on the velocity of the fiycr plate upon impact, 
be it accelerated by a gas gun, magnetic pressure, laser 
light, or explosives. The first-shock end states follow the 
Hugoniot- relation 




Figure 1: (Color online) Theoretical hydrogen Hugoniot curves 
{solid, black: H-REOS.2, short- dashed: H-SCvH-i, long-dashed: 
Sesame-5251) and experimental shock data {grey filled sguares: SNL 
Z-pinch, open squares: modified omega laser, circles: explosives, di- 
amonds: gas gun). The dot- dot- dashed orange curve (S) shows part 
of the Saturn adiabat. 



where initial pressure Pq{pq,Tq) and initial internal energy 
uo{po,To) are derived from th e EOS. Figure [D s hows ex- 
perimental gas gun d ata from Nellis et al. (|l983l ). Sandia 
Z machine data from lKnudson et al. ( 20041) . the modifie d 
Omega laser data from Knudson and DesiarlaisI (|2009[ ). 
and spherical compr ession data using explosives from 
Boriskov et al.l (|2005l ). By scaling the initial density, the- 
oretical hydrogen EOSs can reasonably be compared to 
deuterium experimental data, although some differences 
between D and H become non-negligible in the molecu- 
lar region, which are probed by the gas gun data, due 
to differences in the molecular vibrational states; see 



- uo ^ 0.5{P - Po){po' - p-^) , 



(6) 



Hoist et al.l (|2012|) for a detailed discussion. Figure [T] also 
shows th e theoretica l hydr ogen Hugoniot curves for H- 
RE0S.2 ( Hoist et all . I2OI2I with additional data points 
by A. Becker pers. comm), for the H-SCvH-i EOS, and for 
the H-Sesame EOS. 

Obviously, the theoretical hydrogen EOSs differ sub- 
stantially from each other. Along the Hugoniot states, the 
Sesame EOS, which precedes even the gas gun data, is rel- 
atively stiff at P < 1 Mbar but the most compressible one 
at '^ 5 Mbars. Conversely, the ab initio H EOS is relatively 
compressible below 1 Mbar compared to the Sesame and 
SCvH-i EOS, and to the spherical-shock compression data, 
with a maximum compressibility at P ~ 0.5 Mbar where 
dissociation occurs. At higher pressures of 1 to 3 Mbar, 
the ab initio EOS runs nearly through the experimental 
central values which indicate a low compression ratio of 



4.25, whereas SCvH-i EOS agrees well with the data up to 
0.8 Mbars but then turns to a large maximum compress- 
ibility at ~ 1 Mbar where in the underlying SCvH-ppt 
EOS the plasma phase transition occurs. These proper- 
ties lead to systematic differences in the resulting Saturn 
models. 



3. Results 

3.1. Structure models with different equations of state 

In this Section we focus on the effect of the different 
EOSs (Sesame, SCvH, LM-REOS) on the resulting Saturn 
models. The presented models have been calculated for 
Yi = 0.18, 2n/ijj =10h 39m, and the default values of the 
other observational constraints as given in Table [T] Figure 
[2] shows the results. 

For each of the considered mantle EOS, the parameters 
Zi, ^2, Afcoroj and Mz behave similarly with Pi_2- The 
deeper the layer boundary, the higher become Zi and Z2 
and, as a response, the lower becomes Mcore, while Mz 
remains nearly constant within ~ SM®. Zero-core mass 
models are possible simply by putting the layer boundary 
sufficiently deep into the planet. The maximum core mass 
is determined by the condition that neither Zi nor Z2 must 
become negative. It is 15Mg for rocky cores and about 
2OM0 for water cores. 

For given Pi_2, water cores are up to 50% heavier than 
rock cores because low-density material requires a larger 
volume. Displaced envelope material must be replaced by 
core material, hence the core becomes heavier. Note that 
water core models that approach the limit Afcorc — ^ have 
been computed for LM-REOS mantle EOS only, but they 
exist for the other mantle EOS as well. 

The core mass responds weakly to the surface tempera- 
ture. The found slight decrease by < 2M0 can intuitively 
be explained by the colder and thus denser envelopes. 
However, the colder outer envelope requires AZi > 0.01 
less heavy elements to match J4, but up to AZ2 — 0.05 
more heavy elements in the inner envelope. Thus the 
found insensitivity of A/coro, and also of A/^, to the surface 
temperature appears to be a more complex compensation 
of different effects in Saturn. 

The general behavior of the solutions can be under- 
stood with the help of Figure [3] It shows the sensitiv- 
ity of the gravitational moments to the internal mass dis- 
tribution as parameterized by the contribution functions 
C2m see Eq. ^. While they have been computed for a 
particular model, i.e. the LM-REOS based Saturn model 
S12-3a as highlighted in Figure [2j:, their properties are 
the sam e for all three-layer m odels. As it is well knowm 
(jZharkov and TrubitsvnI . 119781 ) the higher the order of the 
gravitational harmonic, the farther out are the locations of 
mean and maximum sensitivity and the more pronounced 
is the latter one. J2 and J4 are most sensitive at pressures 
of ~ 0.5 and 0.1 Mbar, respectively. At the 1 Mbar level, 
the sensitivity of J4 has dropped to ~ 30% of its maximum 
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Figure 2: (Color online) Resulting Saturn structure models for dif- 
ferent equations of state for the mantle {solid, black: LM-REOS; 
long-dashed, blue: Sesame EOS; dashed, green: SCvH-i EOS) and 
for the core {thick lines: rocky cores, symbols connected by thin 
lines: water cores), and different surface temperatures {open sym- 
bols: Ti bar = 135 K). For given other input parameters (see Section 
l2.1l for details), each displayed model is uniquely defined by its P1-2- 
value (x-axis). Note that the two upper panels have a common j/-axis 
which changes scale a,t Zi = Z2 = 7 X Zg^i^,., where Z^^i^^ = 1.5%. 
One model is highlighted for each EOS by a filled circle in panel (c). 
The red diamond shows the measured 9.12X solar C/H ratio. 



50 



15 10 



P (Mbai-) 

5 3 2 



1 0.5 



0.1 0.01 



^6 

u 

c 
_o 

o 

c 

54 

a 

C 2 

o 

o 



1 


1 III II 


1 1 


- 


mean values — > O O 


o o 


- 






- 




//C9 1 - 
/ "il 


- 






■,/] 


/ / • 

Rcore ^^^''^ ^^^r/^--'~~^ 


i 1 X 111 





0^-=^— ' 



0.2 0.4 0.6 

r(Rm) 



0.8 



Figure 3: (Color online) Contribution functions of the gravitational 
harmonics J2 (thick solid, blue), J4 (short-dashed, red), Jg (dot- 
dashed, green), and Jg (long-dashed, orange) for the Saturn model 
S12-3a. The bottom x-axis scales linearly with mean radius; the 
top X-axis shows the radii where the pressures of 0.01-50 Mbar oc- 
cur. Layer boundaries are clearly seen in the function cq (thin solid, 
black). Diamonds show the radius where half of the final J2n value 
is reached. 



value. J4 is almost insensitive to the mass below 3 Mbars, 
a typical pressure for the outer/inner envelope boundary 
of LM-REOS based Saturn models. Therefore, Zi changes 
little with Pi_2 for P1-2 > 1 Mbar. As the sensitivity of J2 
is similar to that of J4 and just slightly shifted to higher 
pressures, the computed J2 value is affected by the value 
of Zi as well. But because J2 is to be adjusted by the Z2 
value, i.e. by the mass distribution interior to Pi_2 where 
its sensitivity is low, small changes in J2 for adjustment 
require strong changes in Z2, see Figure l^b- 

LM-REOS based models. Using LM-REOS, we find solu- 
tions for 1.2 < Pi_2 < 5 Mbars. For all of these models, 
Zi is nearly constant no matter what the values of Z2 
and Mcorc are. This is because with the layer boundary 
so deep inside, J4 is little sensitive to the mass distribu- 
tion there. As Z2, in contrast, covers a wide range of 0- 
60%, there are solutions with Zi = Z2, unlike for Jupiter 



( Nettelmann et al.l . 12012 ) 



For LM-REOS based models the layer boundary has to 
be put rather deep inside the planet in order to nearly 
suppress the influence of the matter interior to Pi_2 to the 
values of J2 and J4. Otherwise, say for Pi_2 < 1 Mbar, the 
rise in J2 could no longer be compensated for by a lower 
Z2 value because Z2 already goes to zero. This behavior 
is a direct consequence of the higher compressibility of the 
H EOS at sub- Mbar pressures (Figure [IJ. 

The total mass of heavy elements is 16-20 Mq, less 
than predicted for Jup iter with LM-REOS (-- SOM®, 



than Jupiter's (-^ 6 Zq). 

Sesame based models. Sesame EOS based models have 
the layer boundary between 0.2 and 2.5 Mbar. With 
Pi-2 ^^ 0.2 Mbar, the sensitivity of J4 rises strongly (Fig- 
ure 121) so that Zi has to decrease rapidly in order to not 
produce too high | J4I values. There is no overlap of the 
functions Zi{Pi_2) and 2'2(Pi_2), i-C Sesame EOS based 
Saturn models require a higher metallicity in the inner 
than in the outer envelope. Because of the stiffness of the 
H-Sesame EOS up to 1 Mbar (Figure [IJ, the layer bound- 
ary can be farther out and the Zi values as well as the total 
mass of heavy elements can be up to 2.5 times higher than 
for LM-REOS based models. In fact, the layer boundary 
m,ust be farther out because the rise in Z2 for increasing 
Pi_2 is accompanied by a rapid decrease in Mcoro, which 
must not become negative, a direct effect of the higher 
compressibility of the H-Sesame EOS at higher pressures 
up to 10 Mbar (Figure p. 

SCvH-i based models. From the gross compressibility be- 
havior of the Hugoniot up to 100 GPa one would expect 
the SCvH-i EOS based solutions to fall between those for 
the other two EOSs. This in fact happens with respect 
to the parameters Z2, Mcoro, and Mz- In particular, the 
just slightly higher compressibilities of SCvH-i compared 
to Sesame EOS below 100 GPa suggest just slighty lower 
Z\ values than for that EOS. However, the Z\ values are 
relatively higher than for the Sesame EOS based models. 
This can only be due to higher temperatures along the 
SCvH-i Saturn adiabat, as higher temperatures at given 
pressure level reduce the mass density, which allows for 
more heavy elements to be added. We will encounter the 
same argument again in Section W72k In contrast to SG04, 
we do not find models with Z\ = Z2. Presumably, this 
stems from the application of the more recent, accurate 
J2n data. 

Both Sesame and SCvH-i EOS based models have Mz = 
25-30 M^, a similar amount as Jupiter may have (~ 15- 
40Me; SG04). With Zp = 17-21 Zq, this is a larger en- 
richment than in comparable models for Jupiter (3-8 Zq). 



Note that SG04 use an earlier value Z, 



('89) 



0.019. 



Nettelmann et al.l 120121) . However, Saturn's Mz corre- 
sponds to an overall enrichment Zp = 10-15 Zq, higher 



Internal profiles. For each of the EOS, one selected inte- 
rior profile is shown in Figure |4l The temperatures along 
the SCvH-i Saturn adiabat in the outer envelope are in- 
deed higher than for the Sesame adiabat. For the LM- 
REOS based adiabat, the onset of dissociation occurs at 
^ O.TPsat where the temperature gradient fiattens. A typ- 
ical value of the density is 2 g cm~^ in the inner envelope, 
15 g cm~'^ in a rocky core, and 8 g cm~'^ in a water core, 
while in the outer envelope the density changes by four 
orders of magnitude. 

Concluding, for each of the considered EOS, similar val- 
ues of the metallicity in the outer envelope (e.g., 3x solar), 
in the inner envelope (e.g. 10-30 x solar), and of the core 
(e.g., 0-15 AIq) for rocky or 0-18 M® for water cores) are 




0.5 
r (Rsat) 

Figure 4: (Color online) Interior density and temperature profiles of 
the Saturn models that are highlighted in Fig. [2}; by filled circles.; 
solid, black: model with LM-REOS, dashed, green: SCvH-i, and 
long-dashed, blue: with Sesame EOS. These models have an isother- 
mal rock core. Between 0.95 and 1 -Rsat, the density changes by 
three orders of magnitude (not displayed), and between 0.95 -Rsat 
and the boundary to the inner envelope, by one order of magnitude. 
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possible. It depends mainly on the position of an inter- 
nal compositional gradient (for instance in form of a layer 
boundary), which values are adopted. In contrast, the dif- 
ferent EOSs require different locations of that gradient, 
and LM-REOS predicts a lower total mass of heavy ele- 
ments than the other two EOSs. 

3.2. LM-REOS based structure models with different he- 
lium abundances, rotation rates, and J4 values 

The resulting atmospheric heavy element enrichment of 
'^ 3x solar of the LM-REOS based Saturn models is rather 
low compared to the measured carbon enrichment of 9x 
solar. Therefore, we investigate in this Section qualita- 
tively the effect of the assumed atmospheric helium abun- 
dance, of the rotation rate, and of the J4 value on the 
resulting outer envelope metallicity. As we are aiming to 
get higher Zi values than for the models in Section [XTl we 
only consider higher | J4I values, up to the 6(7 observational 
uncertainty (any \J2n\ increases with the mass density in 
the sensitive region), lower Yi values {Zi decreases linearly 
with Yi) of respectively 0.16 and 0.1; but we consider a 
higher rotation rate (to face a possible reality), although 
it will require a reduction in the heavy element content of 
the outer part of a planet. For all of these models, the 
transition pressure is at 300 GPa and the core consists 
of rocks. At temperatures between 140 and 300 K, the 
analytic Van-der-Waals-gas EOS is used for H2 and He, 
and the ideal gas EOS for molecular H2O, but this has no 
relevant effect on the resulting adiabats. 

Figure [S] shows the resulting enrichment factors for the 
metallicity in the outer envelope -to be compared with the 
(potentially) measured abundances- and also in the inner 
envelope. 



Figure 5: Resulting outer {solid lines) and inner envelope metal- 
licities [dashed lines) in terms of the solar metallicity Zq = 0.015 
for different He/(He-|-H) mass ratios of 0.10, 0.16, 0.18 and different 
periods of rotation of lOh 39m and lOh 32m as labeled, and differ- 
ent given J4 values up to the 6(t uncertainty (x-axis). Some of the 
numerically possible models with Z2 close to zero are not displayed. 



Obviously, the influence of J4 on Zi is weak. Only for 
extremely high | J4I values of 9.6 x 10^^ the effect becomes 
of same size as the effect of a slower rotation by 7 min- 
utes, which is AZi - +1Zq. Note that AZi/Zi > 20% 
constitutes a reather big amplification of the 1% relative 
uncertainty in the rotation rate. 

In case a high rotation rate of lOli 32m should proof 
true, the heavy element enrichment would decrease down 
to 2x solar for a helium abundance Yatm = 0.18. Note 
that we do not adjust the planet radius to the rotation 
rate. Interestingly, the rotation rate mainly affects the 
inner envelope metallicity. This can be explained by the 
response of J2 on a change in Zi and Z2 as described in 
[Appendix B[ 

The biggest effect on Zi can be achieved by lowering 
the helium abundance, where AZi ^ 1.5% of heavy el- 
ements can be added for AYi ~ 2% of helium. It is 
|AZi|/|AYi| < 1 because the less helium there is, the 
higher the specific heat of the material, and thus the lower 
must be the temperature along the adiabat to keep the 
entropy constant. Colder adiabats give denser envelopes 
and thus allow for less heavy elements to be put into. 

Despite the wide range of considered parameter values, 
with .^1^6 Zq the possible atmospheric heavy element 
abundance remains clearly below the 9fold enrichment of 
carbon. Our LM-REOS based models therefore predict 
0/H to be less than 9x solar. 
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Figure 6: (Color online) (a) Nondimensional axial moment of 
inertia — core mass relation for a subset of the structure models of 
Fig. [2] using the different EOSs LM-REOS {solid, black), SCvH-i 
{dashed, green), Sesame {long-dashed, blue) with rock cores {lines) 
or water cores {s ymbols); (b)^ : Core mass range and A-range; thin 
dotted: from Ref. JHelledl bOllh for a period of rotation of 10h39m24s 
and scaling by R^q (R. Helled, pers. comra. 2013); thick dotted: 
same as thin dotted but sealed by iJmcan; Arrows (vertical position 
has no meaning) Radau-Darwin approximation, see text for details. 



3. 3. Moment of Inertia 

Helledl ( 20111) showed that an observational determina- 



tion of Saturn's axial moment of inertia, C would im- 
pose an additional eonstraint on Saturn's core mass, and 
that the neeessary measurements can be provided by the 
Cassini extended extended mission. Figure [5] presents 
our results for Saturn's nondimensional moment of inertia, 
A = C/MsatRm^ for a representative subset of the models 
of Fig . [2l As pointed out by and in agreement with lHelled 
( 20111 ). we find different A values for different core mass 
values for models that all meet the observed J2 and J4 
values within their tight la bounds. In particular, for a 
fixed core EOS (rocks or water) and mantle EOS (SCvH, 
Sesame, or LM-REOS), and not too small core mass val- 
ues (> IOM0), the relation between Afcoro and A becomes 
unique. However, due to the uncertainties in the core and 
envelope EOS, a measurement of A could further, but not 
unambiguously constrain the core mass. Moreover, our 
physical EOS based, adiabatic models yield a very narrow 
possible rang e for A, that e ven does not overlap with the 
prediction bv iHelledl (2011), see Figure |B}d. Therefore, a 
measured A value would definitely be of great value for 
discriminating between competing Saturn models. 

As moment of inertia measurements for gas giant plan- 
ets are challenging, the Radau-Darwin approximation is 
often used for an estimate. It expresses A in terms 
of the easier accessible J2 and the expansion parame- 
ter m = uj^Rm /GM, and becomes exact only in the 
limit of constant density bodies. Using i?ni = 58201 km, 
P = 10h39m, and the values of Table [TJ we calcu- 



late A(^^) = 0.2283 (solid arrow in Figure [6b)- The 

I — I J k 

more consistent values of Table 1 in IHelledl (|2011l l sug- 
gest A^^^) = 0.2292 (dashed arrow in Figure ib). Inter- 
estingly, our interior model based values are ~ 3% higher 
than A^^^^ indicating that higher-order deformations ( J4) 
play a non-negligible role in the internal mass distribution. 
Notqj that different scalings of A by either the equatorial 
radius as in IHelledl (|201ll) (R. Helled, pers.comm. 2013), 
or by the mean radius changes the value of Saturn's A by 
7.2%. 

3.4. Homogeneous Evolution 

Homogeneous evolution implies a constant mean molec- 
ular weight in every mass shell of the planet with time, 
while inhomogeneous evolution also allows for vertical 
mass transport such as He rain or core erosion. In the 
considered case of homogeneous evolution, we get cooling 
times rsat of 2.56 Gyr for the LM-REOS model, 2.36 Gyr 
(SCvH-i EOS model), and 2.31 Gyr (Sesame EOS model). 

Neglecting the three contributions from angular mo- 
mentum conservation, change of the energy of rotation, 
and from the time-dependence of the irradiation yields 
0.05 Gyr longer cooling times for Saturn. While with 
TSat ^ 2.5 Gyr, the cooling time comes out significantly 
shorter than the age of the solar system of 4.56 Gyr — 
commonly believed to also be the age of the planets within 
an uncertainty of a few Myr accord ing to circumstellar 



disk observations (IStrom et al.l . Il993f ) — it is in agreement 



with pr evious calcula t ions. For instance, using the SCvH- 
i EOS. iGuillot et al.l (Il995h compute 2.6 ± 0.2 Gyr and 
Fortnev and HubbardI (J2003r ) 2.1 Gyr for an adiabatic, ho- 
mogeneously evolving Saturn. 

Our results show once more that the short cooling 
time of a homogeneously evolving Saturn is essentially 
independent on details of the model assumptions such 
as the size of the core. In other words, we con- 
firm the well-known evidence for a real excess luminos- 
ity compared to the predicted luminosity from homoge- 
neous evolution (|Pollack et al.l . Il977l : ISaumon et al.Lll992 : 
Fortnev and HubbardI . l2003l ). 



4. Summary and Discussion 

^.1. Gross features 

We have applied different EOSs to compute structure 
and evolution models for Saturn within the standard ap- 
proach of a layered interior with only few layers that cool 
down homogeneously with time. Because of the large ap- 
plied input parameter space we have selected combina- 
tions of parameters that we believe yield a reliable esti- 
mate of the overall uncertainty in key internal structure 
parameters. In particular, with LM-REOS and assuming 
Fi = 0.18 and 2tt/lo =10h 39m we find Mcoro = 0-17Me, 



^modified after acceptance of this paper 



Mz = 16-20Me, Zatm < 3x solar, and ts =2.6 Gyr; with 
Sesame-EOS we find Mcore = 0-20Me, Zatm < 7x so- 
lar, Mz = 26-30Me, and ts =2.3 Gyr, while SCvH EOS 
based models have values in between. The value of Zatm of 
the LM-REOS based models can be lifted up to a factor of 
two if Yi is lowered down to 0.10. With ts = 2.3-2.6 Gyr, 
the cooling time is significantly shorter than the age of 
the solar system, pointing to a failure of the cooling model 
that is to be sought beyond the uncertainties in Saturn's 
composition or the EOS. 

We encourage measurements of Saturn's moment of in- 
ertia and of the atmospheric abundance of helium and oxy- 
gen for discriminating between the wide range of possible 
Saturn models and for probing the underlying EOS in cer- 
tain pressure ranges. 

4.2. 0:H ratio 

The 0:H ratio could in principle be derived from bright- 
ness temperature measurements at ^ 1 m wavelengths us- 
ing LOFAR (D. Gautier, pers. comm.), which is a set of 
ground-based radio antennas in western Europe. From a 
measured 0:H, in addition to the already measured C:H, 
the atmospheric mctallicity ^atm can be estimated and 
compared with the Zi values of the Saturn structure mod- 
els. For this purpose, we show possible 0:H-Zatm relations 
in Fig. [71 For a given 0:H, Zatm has been computed as the 
sum of the heavy element particle abundances times their 
atomic weight, divided by the sum over all element abun- 
dances times their weight, where we assume He:H = 0.052 
(>^atm = 0.18), C:H= 0.912X solar, and 2x (or 4x, 6x) so- 
lar abundances of the elements {N, P, S, Ne, Ar, Kr, Xe, 
Mg, A l, Ga}, using the solar system abundances of lLoddersI 
(|2003l ). As Jupiter's atmosphere is observed to be strongly 
depleted in Ne, which may a lso be the case for Saturn if 
caused by He sedimentation ( Wilson and Militzerl . |2010| ). 
we also compute Zatm with Ne:H=0. 

Obviously, models with Zi =~ 3Zq would imply a low 
0:H of only 2x solar; Zi ^ 6 x Zq, the upper limit of the 
LM-REOS based Saturn models, would imply 0:H = 6-8 x 
solar. Higher values are possible with the Sesame or the 
SGvH-i EOS. Goncluding, a measured 0:H could tremen- 
dously help to discriminate between the various Saturn 
models. 

4.3. Core mass and formation 

The possible core mass values of Jupiter and Saturn per- 
sistently attract attention, as the hope to infer the planet 
formation process from the " face value" of the present core 
mass continues to exist. Leading candidates for possible 
formation processes are the core accretion scenario, where 
the g aseous envelope is accreted onto a heavy-e l ement 
core (lAlibert et al. . 20051 : lOodson-Robinson et al.L l20ld : 



2012li 



Kobavashi et al.L 120121 ). and the gravitational disk insta- 
bility scenario, where the planet wo uld form through self- 
contr action of a gaseous cloud (e.g. iHelled and Schubert , 
20081) . 
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Figure 7: Relation between the atniospcric 0:H ratio and the at- 
mosheric metallicity, Zatm in solar units, for which different element 
abundances are assumed: elements apart from C and O are 2 times 
solar {thick solid), 4x solar (dashed), and 6x solar (dot- dot- dashed). 
Thin lines: the same, respectively, but with zero-Neon abundances. 
C:H is 9.12x solar and 0:H is displayed on the x-axis. 



Let us assume the initial core mass of Saturn was 
> 15M0. What could that teUs us? Such a 
heavy cor e exceeds the maximum c ore mass of 8Mq 
found by iHelled and Schuberti (|2008[ ) and thus would 
rule out a disk-instability-kind formation for Saturn. 
On the other hand, it would allow for a comfort- 
ably short time scale for core accretion formation 



( Dodson- Robinson et al.L l2010t iKobavashi et al.l 12012 ) 



even for a non- zero abundance of grains in the protoplan- 



etary envelope (iDodson-Robinson et al 



201(1) . 

If the initial core mass of Saturn was ~ 7Mq, then both 
formation processes could have let to the present Saturn. 
Given Saturn's high total heavy element mass of 16~3OM0, 
disk-insta bility formation would indicate a massive proto- 
solar disk (JHelled and Schuberti . 120091) , as well as an early 
presence of planetesimals be fore formation was completed 
(JHelled and Schuberti . liooj) . 

In the case of an initial core mass of < IMq, again both 
scenarios could be possible. Here, core accretion formation 
would require the absence of grains in the envelope in order 
to reduce the gas press ure induced by warm t emperatures 
in an opaque medium ( Hori and Ikomal l2010f ) . 

As main paths toward better constrained core mass 
properties, we support Helled's (2011) suggestion of a mo- 
ment of inertia measurement, and emphasize the need for 
a better understanding of Saturn's envelope structure. 

4.4- Three-layer models and helium rain in Saturn 

The over and over confirmed finding of a too short cool- 
ing time for homogeneously evolving Saturn models, re- 
gardless of details in the structure models, or of the under- 
lying EOS used, and the simultaneously repeated mention- 
ing of a high likelihood for an inhomogene ously evolving 
Saturn as a result of He sedimentation fe.g.. |Pollack et al 



10 



1977t IStevenson and Salpeteij. ll977L ISaumon et all Il992 



1995t iFortnev and Hubbardl . |200 



pt all 
20031) 



Guillot et all . 119951: IFortnev and Hubbard I2003D cast 
doubt on the usefulness of Saturn models that ignore this 
process. 

The most recent theories of H/He phase separation pre- 
dict demixing under relevant planetary conditions to oc- 
cur when hydrogen metallizes under high pressures while 



the helium atoms are still neutral ( Morales et all l2009l : 



Lorenzen et all l2009i 12011 ). In Jupiter, metallization oc- 
curs rather far out in the planet at ^ . 9 Rj where the 



pressure is only 0.5 Mbar (JFrench et all 120121) . but the 



lowest pressure that could be ac hieved in H/He demix- 
ing calculations so far is 1 Mbar (JLorenzen et al.l . 12011 ). 
Indeed, the ab initio data based H/He phase diagram 
suggests demixing in Saturn at least within 1~5 Mbar 
( Lorenzen et all 12011 ). which corresponds to the region 
beneath ~ 2/3 i?sat- Let us thus divide Saturn's man- 
tle into three regions: an outer region down to 1 Mbar 
(^ 2/3i?sat), an innermost helium-rich region where the 
sedimented helium di ssolves again in its surroundi ng, pos- 
sibly a helium-layer ( Fortnev and Hubbardl . 120031 ). and a 
middle region where the P—T conditions favor H/He phase 
separation. 

In the middle region, the helium abundance at given 
P, T will reduce until the remaining He atoms become 
miscib le: the helium abundance will follow the demixing 
curve (jStevenson and Salpeteii . Il977l ). In P-p-space, such 
an inhomogeneous region could simply be described by a 
smoothed layer boundary. In P-T-space, the effect of an 
inhomogcneity could be significant, as it may be corre- 
lated wnth_a_stron£lY_su£eradf£i^ temperature gradi- 
ent (jStevenson and Salpeteii . 119771 ). which would require 
higher envelope metallicities and affect the derived core 
mass (JLeconte and Chabrieil l2012n . Also important may 
be the induced differential rotation fro m sinking d r oplets 
that conserve their angular momentum. ICao et al.l (|2012l) 
showed that a tiny differential rotation would suffice to 
drive the magnetic field and explain its observed dipolar- 
ity. 

In case the sedimented helium dissolves in the lower re- 
gion, that part is well represented by an inner envelope 
with enhanced abundance as in our models. But in case 
the helium rains down to the core, a preferred solution to 
expla in Saturn's high luminosity ( Fortnev and Hubbardl . 
20031) , the upper part of what we count to be core mate- 
rial may in fact be helium. Thus, Saturn's maximum core 
mass could be lower than predicted by our models. 

In the outer region we would mainly see the depletion 
in helium, because helium atoms from the upper regions 
are transported repeatedly by convection down into the 
immiscibility region, where a fraction of them gets lost 
into the deep through sedimentation. Whatever happens 
to the compositional and temperature gradient deep in 
the planet, the outer, miscible envelope should remain ho- 
mogeneous and adiabatic. Therefore, our results for the 
atmospheric helium and heavy element abundances are 
certainly caused or influenced by He-sedimentation, but 



are not expected to alter when He sedimentation would be 
explicitly accounted for. 

^.5. Why would we expect a (dis) continuous heavy element 
distribution? 

Standard three-layer models with heavy element discon- 
tinuity benefit from two additional parameters (Z2, P1-2) 
that can be used to fit the observational constraints. How- 
ever, those models would even more benefit from a phys- 
ical justification. At least, the assumption of a sharp 
layer boundary between two convective, homogeneous lay- 
ers offers the advantage of a self-consistent picture, where 
upward particle transport across would be inhibited but 
not necessarily the heat transport. In contrast, a contin- 
uous, inho mogeneous heavy element di stribution as sug- 
gested by iLeconte and Chabrierl (|2012i ) might lead to a 
semi-convective boundary layer with reduced heat trans- 
port. Whether such a picture can explain the luminosities 
of both Jupiter and Saturn remains to be shown. 

A possible origin for an inhomogeneous heavy element 
distribution could be the erosion of an initially big core, 
with subsequent small-sca le layer formation wit h sma ll 
compositional gradients (JLeconte and Chabrieil I2OI2I ). 
and finally merging of the layers to a single one as seen in 
hydrodynamic simulations of fluids with both tem perature 
and compositional gradients (jWood et alll2012l ). Thus, a 
high, homogeneous metallicity in the inner envelope could 
be the result of an eroded core where the core material 
went through stages of layer merging until the composi- 
tional gradient got large enough to stop merging with what 
we now see as an outer envelope. 

4-. 6. Summarized main findings 

• The H/He EOS strongly influences the atmospheric 
metallicity, Zatm and the possible position of an in- 
ternal layer boundary, but has little influence on the 
core mass and the cooling time. We find Mcorc = 0- 
20Me and Tgat - 2.5 Gyr. 

• The total mass of heavy elements in Saturn can be less 
(LM-REOS) or equal (Sesame, SCvH-i EOS) to that 
in Jupiter, while the averaged enrichment is larger 
than that of comparable Jupiter models. 

• Our LM-REOS based Saturn models predict 2'atm ~ 
3x solar (^ 6x) for an atmospheric helium abundance 
Yatra — 18% (Yatm = 10%) by mass. The correspond- 
ing predicted maximum 0:H ratio is ~ 2x solar (8x). 

• For Saturn, we calculate a non-dimensional axial mo- 
ment of inertia A = 0.235 to 0.236. 

We thank Daniel Gauticr for illuminating conversations 
on abundance measurements, Johannes Wicht and Ravit 
Helled for interesting discussions, Andreas Becker for com- 
puting the high-temperature extension of the H-RE0S.2 
Hugoniot curve, and Winfried Lorenzen for discussions on 
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SO that Api '^ p^AZi and Ap2 ^ P2^^2- Because the 
value of J2 is an observational constraint, AJ2 must be 
zero under the perturbations Api and Ap2- With Equa- 
tions (|B.1[) and (jB.2p we thus have 

= AJ2 « VLiAZiplr\pi) + VL2AZ2plr\p2) , (B.3) 



hence 



AZ2 



-AZ, (YlI 

\Vl2 



riP2 



(B.4) 



For a typical three-layer Saturn model as shown in Fig- 
ure H r{P,_,) - 0.5i?sat, leading to Vli/Vl2 ~ l/O.S^ = 
23, r{pi) = 0.7i?sat, r{p2) = 0.35i?sat, leading to 
(r(pi)/r(p2))^ - 2\ pi - 0.7gcm-3, p, ^ 2.05 gcm'^ 
leading to {pi/p2f ~ (l/3)^ and thus AZ2/AZ1 « 
—2^/3^ K, —16, in reasonable agreement with AZ2I AZ\ « 
— 10 as seen in Figure [5] 



Appendix A. Erratum Jupiter-paper 



In the Jupiter-II paper by iNettelmann et al.l (|2012i ) on 
Jupiter structure and homogeneous evolution models it 
was stated that including the energy of rotation in the 
thermal evolution leads to a reduced luminosity with time 
and thus a longer cooling time. The latter statement must 
be corrected for a shorter cooling time: the energy of ro- 
tation will not be released by radiation at a later time. 
Therefore, Jupiter's cooling time is decreased by 0.2 Gyr 
(and not prolonged by 0.2 Gyr as stated in that paper). 
The given final value for cooling time of 4.41 Gyr, when in 
addition the time-dependence of Lq is included, remains 
valid. 



Appendix B. Estimate AZ2 as a function of AZi 



We here derive an estimate for the necessary change in 
Z2 in response to a change AZi when J2 is to be kept 
unchanged. According to Equation [1] wc can approximate 
J2 by 

J2 « Vli pi r^ipi) + Vl2 P2 r\p2), (B.l) 

where pi, i = 1,2, is the mean density of layer No. i, and 
Vn its volume. The contribution from the small, central 
core is neglected. According to the additive volume law 
for mixtures we can write p^^ = (1 — Z)p^ ^^ + Zp^ and 
thus 
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